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Abstract. We study spherical Whittaker functions on a metaplectic cover of 
GL(r + 1) over a nonarchimedean local field using lattice models from statistical 
mechanics. An explicit description of this Whittaker function was given in terms 
of Gelfand-Tsetlin patterns in [SJ[T7], and we translate this description into an 
expression of the values of the Whittaker function as partition functions of a six- 
vertex model. Properties of the Whittaker function may then be expressed in terms 
of the commutativity of row transfer matrices potentially amenable to proof using 
the Yang-Baxter equation. We give two examples of this: first, the equivalence of 
two different Gelfand-Tsetlin definitions, and second, the effect of the Weyl group 
action on the Langlands parameters. The second example is closely connected 
with another construction of the metaplectic Whittaker function by averaging over 
a Weyl group action [9|fl0] . 



1. Introduction 

The study of spherical Whittaker functions of reductive groups over local fields is 
of fundamental importance in number theory and representation theory. Recently, in 
two separate series of papers, the authors and their collaborators have studied Whit- 
taker functions on metaplectic covers of such groups. The goal of this paper is to 
introduce a new method for describing such p-adic metaplectic Whittaker functions: 
two-dimensional lattice models of statistical mechanics. In such a model, one defines 
the partition function to be a weighted sum over states of the model. We show 
that there exists a choice of weights for which the partition functions are metaplec- 
tic Whittaker functions. Baxter j2] developed important techniques for evaluating 
the partition functions of lattice models including the so-called "commutativity of 
transfer matrices" and the use of the Yang-Baxter equation. We discuss how these 
methods relate to our descriptions of Whittaker functions and to prior work. 

Two different explicit formulas have been given in [S] and [T7] for the spherical 
Whittaker function on a metaplectic cover of GL(r + l) over a non-archimedean local 
field. The first of these is expressed in terms of a Weyl group action described in [10J, 
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the second in terms of a function on Gelfand-Tsetlin patterns initially introduced 
in [7]. In fact, this latter representation belongs to a family of explicit formulas, one 
for each reduced expression of the long element of the Weyl group as a product of 
simple reflections. Two such reduced expressions in type A are particularly nice, 
and lead to representations of the Whittaker function as sums over Gelfand-Tsetlin 
patterns. In keeping with earlier works, we refer to the two different descriptions as 
"Gamma" and "Delta" rules. The main result of |Bj is a combinatorial proof that 
these two definitions are in fact equal. This equality allows one to prove the analytic 
properties of an associated global object (a multiple Dirichlet series) by applying 
Bochner's convexity principle. 

In the following section we demonstrate that the Gelfand-Tsetlin patterns we are 
concerned with are in bijection with admissible states of the six-vertex model having 
certain fixed boundary conditions. After recalling the description of the metaplectic 
Whittaker function as a function on Gelfand-Tsetlin patterns in Section 3, we use the 
bijection with the lattice model in Section 4 to express both the Gamma and Delta 
descriptions of the Whittaker function as partition functions for certain respective 
choices of Boltzmann weights. 

In Section 5, we take the connection with statistical models further. We show 
that the necessary result for demonstrating the equivalence of the Gamma and Delta 
descriptions may be reformulated in terms of the commutativity of transfer matrices. 
Baxter [2j advocated the use of the Yang-Baxter equation for demonstrating this 
commutativity. In Section 6, we explain how this is carried out in the context of the 
six-vertex model and we speculate about the possibility of such an equation in the 
metaplectic case. 

Finally we discuss the Weyl group action on metaplectic Whittaker functions, 
initially established by Kazhdan and Patterson [13J, which plays a critical role in 
the explicit formulas of [9]. When the degree of the cover is 1, i.e. the linear case, 
the p-adic spherical Whittaker function is essentially a Schur polynomial by results 
going back to Shintani [19]. The Weyl group action is thus closely related to the 
standard permutation action on polynomials in r + 1 variables. In [3], this Whittaker 
function (or equivalently, the Schur polynomial multiplied by a g-deformation of the 
Weyl denominator) is realized as a partition function on a six-vertex model and its 
properties are studied via instances of the Yang-Baxter equation. On the other hand, 
as soon as the degree of the cover is greater than 1, the action looks rather different 
(cf. (I3~l~]) - (l33]) ). Nevertheless, we may ask whether these functional equations may 
also be phrased in terms of transfer matrices and a Yang-Baxter equation, and in 
this final section, we present evidence towards an affirmative answer. 

This work was partially supported by the following grants: NSF grants DMS- 
0844185 (Brubaker), DMS-0652817 and DMS-1001079 (Bump), DMS-0847586 (Chinta), 
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grant H98230-10-1-0183 (Friedberg). 



2. Six Vertex Model and Gelfand-Tsetlin Patterns 



In this section, we demonstrate a bijection between strict Gelfand-Tsetlin patterns 
and admissible states of the six vertex model (or "square ice") on a finite square 
lattice with certain fixed boundary conditions. The boundary conditions on ice were 
known to Hamel and King, who presented bijections between ice and patterns related 
to the symplectic group in [T2]. A treatment tailored to the aims of the present paper 
was given in [I], whose terminology we now recall. 

A Gelfand-Tsetlin pattern of rank r is a triangular array of integers 



a o,o 



a o,i 



Qo,r-l 



a 0,; 



in which the rows interleave: 



^ ay ^ Cli 



The set of all Gelfand- 



Tsetlin patterns with fixed top row is in bijection with basis vectors of a corre- 
sponding highest weight representation of GL(r + 1,C). Indeed, any given top row 
(a ,o, Oo,i5 • • • > a o,r) is a partition which may be regarded as a dominant weight of 
the GL(r + 1, C) weight lattice. Each successive row of a pattern then records a 
branching rule down to a highest weight representation on a subgroup of rank one 
less. We will focus mainly on the set of strict Gelfand-Tsetlin patterns, whose entries 
in horizontal rows are strictly decreasing. In terms of representation theory, these 
patterns result from branching through strictly dominant highest weights. Top rows 
of strict Gelfand-Tsetlin patterns are then indexed by strictly dominant weights A + p 
where A is a dominant weight and p is the Weyl vector (r, r — 1, . . . , 0). 

Now we come to lattice models. The six-vertex model consists of labelings of 
edges in a square grid where each vertex has adjacent edges in one of six admissible 
configurations. This model is sometimes referred to as "square ice" where each vertex 
of the grid represents an oxygen atom and the 6 admissible ways of labeling adjacent 
edges correspond to the number of ways in which two of the four edges include a 
nearby hydrogen atom. If we represent adjacent hydrogen atoms by incoming arrows, 
and locations where there is no adjacent hydrogen atom by outgoing arrows, the six 
admissible states are as follows. 
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We will use a representation consisting of a lattice whose edges are labeled with 
signs + or — , called spins. To relate this to the previous description, interpret a 
right-pointing or down-pointing arrow as a +, and a left-pointing or up-pointing 
arrow as — . We then find the following six configurations. (The index i in the table 
indicates the row to which the vertex belongs, and will be used in later sections.) 



e 



e 




e- 



-e 



e 







0- 

i 





The rectangular lattices we consider will be finite, with boundary conditions chosen 
so that the admissible configurations are in bijection with strict Gelfand-Tsetlin 
patterns with fixed rank r and top row A + p as above. Here A = (A r , . . . , Ai, Ao) 
with Xj > Xj-x for all j, and we suppose that Ao = 0. 

Boundary Conditions. The rectangular grid is to have A r + r + 1 columns (labeled 
through X r + r increasing from right to left) and r + 1 rows. Then with X + p = 
(X r + r, A r _i + r — 1, . . . , 0), we place a — spin at the top of each column whose label 
is one of the distinct parts of A + p, i.e. at columns labeled Xj + j for ^ j ^ r. We 
place a + spin at the top of each of the remaining columns. Furthermore, we place 
a + spin at the bottom of every column and on the left-hand side of each row and a 
— spin on the right-hand side of each row. 

For example, put r = 2, and take A = (3, 2, 0), so that A + p = (5, 3, 0). Then we 
have the following boundary conditions for the ice: 

5 4 3 2 1 

e ffi e © © e 



(2) 



Q Q Q Q Q Q 



© 







©©©©©© 

The column labels are written above each column, and row labels have been placed 
next to each vertex. These row labels will be used in Section HJ but need not concern 
us now. The edge spins have been placed inside circles located along the boundary. 
The remaining open circles indicate interior spins not determined by our boundary 
conditions, though any filling of the grid must use only the six admissible configura- 
tions in the above table. Such an admissible filling of the finite lattice having above 
boundary conditions will be referred to as a state of ice. 
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Proposition 1. Given a fixed rank r and a dominant weight A = (A r , . . . , Ai, 0), 
there is a bisection between strict Gelfand-Tsetlin patterns with top row A + p and 
admissible states of ice having boundary conditions determined by A as above. 

Proof. We begin with a strict Gelfand-Tsetlin pattern. Each row of the Gelfand- 
Tsetlin pattern will correspond to the set of spins located between numbered rows 
of ice, the so-called "vertical spins" since they lie on vertical edges of the grid, as 
follows. To each entry a^j in the Gelfand-Tsetlin pattern, we assign a — to the 
vertical spin between rows labeled r + 2 — i and r + 1 — i in the column labeled ajj. 
(Recall that we are using decreasing row labels from top to bottom as in the example 
02]).) The remaining vertical spins are assigned +. 

It remains to assign horizontal spins, but these are already uniquely determined 
since the left and right edge horizontal spins have been assigned and each admissible 
vertex configuration has an even number of adjacent + spins. We must only verify 
that the resulting configuration uses only the 6 admissible configurations (from the 8 
having an even number of + signs) for the corresponding ice. This is easily implied by 
the interleaving condition on entries in the Gelfand-Tsetlin pattern, which is violated 
if one of the two inadmissible configurations appears. See Lemma 2 of [4] for more 
details. □ 

A simple example illustrates the bijection: 

5 4 3 2 1 



We now discuss the relation between the spherical metaplectic Whittaker function 
on the n-fold cover of of GL(r + 1) over a non-archimedean local field and Gelfand- 
Tsetlin patterns. Such a relationship, described globally, was conjectured in [TJ and 
was established in [5]. Though it is possible to pass from the global result to its 
local analogue (cf. [II]), a direct local proof was given by McNamara [IT], expressing 
a metaplectic spherical Whittaker function as a generating function supported on 
strict Gelfand-Tsetlin patterns. In this section, we recall two formulations of this 
explicit description, following [6] . In Section H] we will explain their translations to 
square ice via the bijection of Proposition [TJ 



(3) 




3. Metaplectic Whittaker Functions and Patterns 
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Let G(F) denote the n-fold metaplectic cover of G(F) = GL(r + 1,F), where F 
denotes a nonarchimedean local field having ring of integers and residue field of 
order q. (There are several related such extensions, but we will use the one in [T3] 
where their parameter c = 0.) We assume 2n divides q — 1, which guarantees that F 

contains the group °f 2nth roots of unity. The group G(F) is a central extension 
of G(F) by yL n : 

l^lin^ G(F) G(F) — ► 1. 
We will identify \i n C F with the group fi n C C of complex ra-th roots of unity by 
some fixed isomorphism. For convenience, we will sometimes denote G(F) as just G, 
and if H is an algebraic subgroup of G, we may denote by H the preimage of H(F) 
in G. 

For details of the construction of the metaplectic group and results about its rep- 
resentations, see [TS] in this volume. Let B(F) be the standard Borel subgroup of 
upper triangular matrices in G(F), and let T(F) be the diagonal maximal torus. 
Then B(F) = T(F)U(F) where U(F) is the unipotent radical of B(F). The meta- 
plectic cover splits over various subgroups of G(F); for us it is relevant that it splits 
over U(F) and over K := G(of), the standard maximal compact subgroup. By 
abuse of notation, we will sometimes denote by K the homomorphic image of K in 
G under this splitting. 

Let s : G(F) — > G be any map such that 7r o s is the identity map on G(F). Then 
the map o : G(F) x G(F) — > \i n such that 

s(gx)s(g 2 ) = a(gi, g 2 )s(g 1 g 2 ) 

is a 2-cocycle defining a class in H 2 (G(F), /x n ). A particular such cocycle was con- 
sidered by Matsumoto [16], Kazhdan and Patterson [13] and Banks, Levy and Sep- 
anski [lj. By these references, such as (13] Section 0.1, the map s may be chosen so 
that the restriction of a to T(F) is given by the formula 

(4) *(( ],( ))=II(^«i)- 

V V U +1 J \ u T+1 J J i<i 

The cocycle a also has the property that a(u,g) = a(g,u) — 1 if u G U(F), and so 
the restriction of s to U(F) is a homomorphism to G. 

We will call a representation tt of G or any subgroup genuine if 7r{(g) = (^{g) when 
£ G fi n . Recall that T denotes the inverse image under 7r of the maximal torus T(F). 
The center Z(T) of T has finite index. Since T is a Heisenberg group, characters 
of Z(T) parametrize the irreducible genuine representations of T as follows. Let \ 
be a quasi-character of Z(T) that is unramified. This means that it is trivial on 
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the maximal compact subgroup of this abelian group. Let A be a maximal abelian 
subgroup of T. Extend x t° A then induce it to T. By Clifford theory, the resulting 
representation i(x) of T is genuine, irreducible, and finite-dimensional. It does not 
depend on either A or the extension of \ to A. 

We extend i(x) to the inverse image B of B(F) in such a way that s(U(F)) acts 
trivially. We then consider the representation of G obtained by normalized induction. 
We call the vector space of the resulting representation I(x)- It has a one- dimensional 
space of -fT-fixed, i.e. spherical, vectors. 

Let : F — > C be an additive character that is trivial on Of but on no larger 
fractional ideal. Given a spherical vector / in I(x), we may construct the spherical 
metaplectic Whittaker function via the integral 



(5) 



W(g) 



f(w s(u)g)ijj(u)du, 



U(F) 



where ip is the character of U(F) given by 



//l Xl,2 ■■■ 
1 X 2 ,3 



\ 



\ 



Xl,n\ \ 
1 II 



and wo is a representative in K of the long Weyl group element. Strictly speaking, 
W(g) as we have defined it is an i(x)~ valued function and should be composed with a 
natural choice of linear functional on i(x) to obtain a complex- valued function. For 
brevity we will sometimes refer to W(g) as the metaplectic Whittaker function. 

Note that x on may be described by a choice of r complex numbers s = 

(si, . . . , s r ). The transformation property W(s(u)gk) = ip(u)W(g) for all u G 
U(F),k G K, implies that it suffices to determine W on the inverse image of the 
torus T(F). Moreover, since W is genuine, it is sufficient to specify W on s(T(F)). 
Given A = J2i ^i^ii where are fundamental weights, let t\ be the element of the 
torus T(F) given by 

/pXi+X'2-i hA r \ 



tx 



p 



A 2 +-+A r 



V 



p 



where p is a uniformizer for Of- 



1/ 



Let = s(t\). Due to our assumption that F contains the 2n-th roots of unity, 
i.P:P)n = 1 an d by (j3J it follows that t A+M = t A t M . 

It is not hard to show that W(t\) = unless A is a dominant weight. 
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Given any dominant weight A, the metaplectic Whittaker function W(t^) may 
thus be expressed as a Dirichlet series in r complex variables s — (si, . . . , s r ) of the 
form 

(6) Z(s;X)= J2 H(p k \...,p kr ;X)q kl{1 - 2si)+ - +kAl - 2Sr \ 

k={k 1 ,...,k r )eN r 

where q denotes the cardinality of the residue field Op/p. We now describe the 
function H(p kl , . . . ,p kr ; A). See [T7] for further details. 

The positive integer n will continue to denote the degree of the metaplectic cover. 
We define the Gauss sum 

g(a,b)= (u,p) b n ip (p a ~ b u) du, 

Jo* 

where (•, -) n denotes the n-th power Hilbert symbol, and we normalize the Haar 
measure so that fi(op) = 1. As a further shorthand, for any positive integer b, we set 

(7) g(b)=g(b-l,b), h{b)=g{b,b). 

Note that for a fixed base field F, these values depend only on b mod n. If n divides 
b, in particular if n — 1, we have 

(8) g(b) = ~, h{b) = i-- q 

We caution the reader that the g-powers that appear in the g{a,b) are normalized 
differently than in the previous works [HH7] ; these are the functions denoted and 
in j6]. The function h is a degenerate Gauss sum whose values may be made explicit, 
while (if n\b) g(b) is a "true" Gauss sum. 

Any strict Gelfand-Tsetlin pattern T with entries indexed as in (JTJ, we associate 
a weighting function 7 to each entry a^j with % > 1 as follows: 

(gipij) if dij = a^i^x, r 

KKi) if a i-iJ ^ a i,j ^ Oi-ij-i, where b itj = ^(oi,f-Qt-i,t)- 

1 if dij = Qi— ij, i=j 

Then we define 



(10) Gr ( T )=nn^ 

i=l j=i 

If T is a Gelfand-Tsetlin pattern that is not strict, we define G r (T) = 0. We also 
define 

r 

(11) k r (%) = (A£(X), . . . , fc r r (T)) where kf (T) = ^ a*,, - a , ; . 
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In particular, note that both G T and fc r are defined using differences of elements 
above and to the right of dij. The superscript T may be regarded as indicator that 
these quantities are defined using such "right-hand" differences. 

We present these definitions in this ad-hoc fashion in order to give a brief and self- 
contained treatment, but in fact they have very natural descriptions when reinter- 
preted as functions on a Kashiwara crystal graph. See [6] for an extensive discussion. 

As an example, consider the Gelfand-Tsetlin pattern ^ in (J3J). Then 
(12) 

(61,1,61,2,62,2) = (1,1, 2) so that G r (%) = h(l)g(2), and (k u k 2 ) = (1, 3). 

Theorem 1. (Brubaker, Bump, and Priedberg [5j; McNamara [1TJ ) Given 
a dominant weight X and a fixed r -tuple of non-negative integers k = (k\, . . . , k r ), 
the function H(p kl , . . . ,p hr ; A) appearing in the p-adic Whittaker function W(t\) is 
given by 

H(p k ;X):=H(p k \...,p k -,X) = ^ G r (T), 

fc r (X)=fc 

where the sum is over all Gelfand-Tsetlin patterns with top row corresponding to X+p 
satisfying the subscripted condition. 

There is a second explicit description of H(p k ; A) in terms of "left-hand" differences 
using functions G A and k A that are analogous to those defined in (flOl) and (TTTT) 
respectively. Assuming that T is strict, set 
(13) 

{g(cij) if dij = di-ij, 3 
h(cij) ifai-ij ^ a itj ^ ai-ij-i, where cy = ^^fa-hi-i ~ a hi) 
1 if dij = cti-ij-i, 1=1 

and define 

r r 

(14) G A (1) = lH[5(a h ,). 

i=l j=i 

If 1 is not strict, define G A (T) = 0. We also set 

i 

(15) fc A (X) = (k A CX),...,k A (<£)) where k A {%) = ^ a , { -i - Oh-i-^+i-j. 

1=1 

The main theorem of [B] is as follows: 

Theorem 2. (Statement A of Brubaker, Bump and Friedberg [6]) Given a 
dominant weight X and a fixed r -tuple of non-negative integers k = (k%, . . . , k r ), 

(16) Yl ° r w= E 

k r (%)=k fc A (X)=fe 
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where the sums each run over all G elf and- Tsetlin patterns with top row corresponding 
to A + p satisfying the subscripted condition. 

As an immediate corollary, we have a second description of the p-adic Whittaker 
function in terms of G A and k A . We refer to these two recipes on the left- and 
right-hand sides of ( |T6l) as the T- and A-rules, respectively. 

In fact, there are many other descriptions for the Whittaker function, though these 
are generally much more difficult to write down as explicitly Indeed, as explained 
in [3l[T5], there exist bases for highest weight representations corresponding to any 
reduced expression for the long element wo of the Weyl group of GL(r + 1) — S r , the 
symmetric group on r letters — as a product of simple reflections ov These make 
use of the Kashiwara crystal graph and are commonly called string bases. Using 
these bases, one may make a correspondence between long words and recipes for the 
Whittaker function (cf. [6J Chapter 2]). From this perspective, the T-rule corresponds 
to the word 

w = eri(er 2 0"i) • ■ ■ (oV^r-l ■ ■■a 1 ), 
whereas the A-rule corresponds to the word 

Wo = ov(ov_i<T r ) • • • (<7icr 2 • • -a r ). 

These two words are as far apart as possible in the lexicographic ordering of all 
reduced decompositions. The proof of Theorem [2] as given in [6] uses a blend of 
combinatorial arguments to give various equivalent forms of the identity ffT6l) as we 
move through the space of long words. We highlight various aspects of the proof in 
more detail now. 

The proof is by induction on the rank r. The inductive hypothesis allows us to 
equate any two recipes for the Whittaker function whose associated long words differ 
by a sequence of relations obtained from a lower rank case. For example, assuming 
the rank 2 case allows us to perform a braid relation o\O20~\ = cr 2 o"iO"2, which could 
be applied to the word corresponding to the T-rule above. After a series of such 
identities, we arrive at two descriptions for the Whittaker function as a weighted sum 
over Gelfand-Tsetlin patterns that agree on the bottom r — 2 rows of the pattern. 
Thus, we may restrict our attention to the top three rows of a rank r pattern. We 
refer to such three-row arrays of interleaving integers, where we fix both the top and 
bottom of the three rows, as "short Gelfand-Tsetlin patterns" and re-index the three 
rows as follows: 

( t t x ■■■ 4_i 4 

(17) t = < d\ 02 CLr-l a r 

I mi vrt2 m r -i 

These two recipes for the Whittaker function will be called G rA (as this recipe uses 
a right-hand rule for the entries and a left-hand rule for the entries rrij) and G AT 
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(where the use of rules is reversed). To be exact, using the definitions in fl9]) and 
(113]) . we have 

r i — 1 r r— 1 

G rA (t)=n 7 (aon*(™i)> and GAr ( t )=n^)n^K0- 

i=i j=i i=i j=i 

Rather than define functions k rA and k Ar on short patterns in analogy to the recipes 
above, it is enough to specify the middle row sum as the other rows are fixed. 

Before stating the reduction, we require one final ingredient. There is a natural 
involution q r on short Gelfand-Tsetlin patterns of rank r, given by acting on middle 
row entries according to 

q r : Oj i — > max(£j_i, m^-i) + min(^, m,) - a, =: a-, 

where if z = we understand that max(4, mo) = £o and if i — r, mm(£ r ,m r ) = 
£ r . This involution q r is used by Berenstein and Kirillov (cf. [H]) to define a 
Schiitzenberger involution on Gelfand-Tsetlin patterns. Brubaker, Bump and Fried- 
berg use the involution q r to give the following reduction of Statement A in [6]. 

Theorem 3. (Brubaker, Bump and Friedberg; Statement B of |6J) Fix an 

(r + 1) -tuple of positive integers £ = (£q, . . . , £ r ), an (r — 1) -tuple of positive integers 
m = (rrii, ■ ■ ■ , m r-i)i and a positive integer k. Then 

£ G rA (t) = G Ar (q r (i)), 

where a\ are the entries of q r (t), k' = ^\ + '^ J - nij — k, and the sums range over all 
short patterns with top row £ and bottom row m satisfying the indicated condition. 

See [8] and Chapter 6 of [6] for a full proof of the reduction from Statement A to 
Statement B. As noted above, the proof of Statement B proceeds through a series of 
additional reductions which occupy thirteen chapters of [Bj. In brief, for "generic" 
short patterns t, the Schiitzenberger involution q r gives a finer equality G rA (t) = 
G Ar (g r (t)), which implies the equality of Statement B summand by summand. By 
"generic" we mean that the entries of the short pattern are in general position - in 
particular, for all i, £i ^ using the notation of (1T7I) . Note that the Schiitzenberger 
involution does not necessarily preserve strictness for all patterns in the remaining 
non-generic cases, and one needs much more subtle arguments to handle these short 
patterns. For such patterns, Statement B is not in fact true summand by summand, 
and one does need to sum over all short patterns with fixed row sum to obtain 
equality. 

As an alternative to establishing Statement B, we mention that one could also 
prove Theorem [2] by computing the Whittaker integral in two ways, mimicking the 
techniques of [5], thus obtaining a proof via decomposition theorems in algebraic 
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groups which respect the metaplectic cover. In subsequent sections of this paper, 
we propose a third way of viewing these theorems using ice-type models, which por- 
tends new connections between number theory/representation theory and statistical 
physics. 

4. Ice and Metaplectic Whittaker Functions 

In statistical mechanics, one attempts to infer global behavior from local interac- 
tions. In the context of lattice models, this means that we attach a Boltzmann weight 
to each vertex in the grid and for each admissible state of the model, we consider 
the product of all Boltzmann weights ranging over all vertices of the grid. Then 
one can attempt to determine the partition function of the lattice model, which is 
simply the sum over all admissible states of the associated weights. In this section, 
we explain how to obtain the metaplectic spherical Whittaker function of Section [3] 
as the partition function of a lattice model with boundary conditions as defined in 
Section [5J 

We make use of the two sets of Boltzmann weights B r and B A . When these weights 
are applied to an admissible state of ice, we refer to the resulting configuration as 
Gamma ice or Delta ice, respectively. In order to indicate which set of weights is 
being used at a particular vertex, we use • for Gamma ice and o for Delta ice. 
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In giving these Boltzmann weights, we have made use of the notation in (J7|). For 
Gamma ice, the constant a equals the number of + signs in the z-th row to the right 
of the vertex •. For Delta ice, the constant a equals the number of — signs in the 
z-th row to the left of the vertex o. In either case, we refer to this constant as the 
"charge" at the vertex. Note by our definitions in (JTj), the Boltzmann weights only 
depend on the charge mod n. The weights B r and _B A also depend on parameters 
Zi, where i indicates the row in which the vertex is found. For Gamma ice, the row 
numbers decrease from r + 1 to 1 as we move from top to bottom as in the example 
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([2]), while for Delta ice, the row numbers increase from 1 to r + 1. These Zi are 
referred to as "spectral parameters." We often suppress the dependence of B v and 
B A on the spectral parameters Zi, 1 < i < r + 1. Let z = (zi, . . . , z r+ \). 

Given an admissible state of Gamma ice (or Delta ice, respectively) S, we define 

(19) Q r (S, z) = J] B r (v), g A (S, z) = J] B A (v), 

where the product (in either case) is taken over all vertices in the state of ice S. 

Proposition 2. Under the bijection of Proposition [IJ with strict pattern X corre- 
sponding to an admissible state of Gamma ice S, then G r (%) as defined in ^^) is 
related to Q r (S,z) in / TJgj) as follows: 

where di(%) is the sum of the entries in the i-th row of the pattern T. 

Similarly, for an admissible state of Delta ice S, G A {%) as defined in [Lfy is related 
to Q A (S,z) in m$ by 



G A (S,z)=G A (Z)z a 1 



doCr)-di(S) di(X)-d 2 (X) 



L (X)-d r (X) d r {%) 



We first illustrate this for Gamma ice with our working example from ([3]) in rank 2. 
The admissible Gamma ice S and its associated Boltzmann weights are pictured 
below. 
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Taking the product over all these weights, we obtain g r (S,z) = h(l)g(2)z^z 2 zf. 



which indeed matches G r ( c £)z 3 
as in ffT2l). 



do (X) -di (X) di (X) -d 2 (X) d 2 (X) 



with T as in © and G r (T) 



Remark 1. The relevant terms in the metaplectic Whittaker function take the form 

as given in Theorem^ However, k r (%) = (k\, . . . , k^) may be easily recovered from 
our fixed choice of highest weight A + p = (£\, - ■ • , £ r , 0) and the row sums dj := cfj(X) 
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used in the monomial above. Indeed, 

k\ = d 1 - (£ 2 + h 4), k T 2 = d 2 - (£ 3 + ■■■ + £ r ),..., k T r = d r . 

Hence, upon performing this simple transformation, we may recover the monomials 
in g 1 " 2 ^ in §20\) from those in Zj appearing in Q r (S, z) of the above proposition. A 
similar set of transformations holds for the Delta rules. 

Proof. Proposition [2] is a consequence of the bijection given in Proposition [TJ We 
sketch the proof for Gamma ice, as the proof for Delta ice is similar. Recall that — 
vertical spins correspond to the entries of the pattern, so the values 7(dij) given in 
(Q should appear in the Boltzmann weights for vertices sitting above a — vertical 
spin. The particular cases of Qj to be used are determined by the vertical spin above 
the vertex in question. We now show that each b it j in (Q matches the charge, the 
number of + signs to the right of the vertex in row i. Equivalently, we must show 
that every spin between column a* j and column aj-ij in row i is assigned a +. So 
suppose that a it j > Oj-ij and let v be the vertex in row i, column a it j, and let w be 
the vertex in row i, column dj-ij. Then the north and south spins for v are (+, — ) 
which, by the six admissible configurations in Gamma ice, forces the east spin to be 
+. All the vertices between v and w have north and south spins (+, +) according 
to our bijection. The east spin + for v becomes the west spin for the neighboring 
vertex v' to the right of v, forcing the east spin of v' to be + as well. This effect 
propogates down the row, forcing all row spins between v and w to be +. Finally we 
must show that the spectral parameters for Q are given by differences of consecutive 
row sums. This is Lemma 3 of [3]. □ 

Given a fixed set of boundary conditions for the vertex model and an assignment 
B of Boltzmann weights associated to each admissible vertex, we refer to the set of 
all admissible states 5asa "system." Given a system (5, its partition function Z(&) 
is defined as 

(21) Z{6):=Z{&,z) = Y,B(S,z), with B(S, z) : = J] B(v), 

see ves 
where this latter product is taken over all vertices v in the state S. In particular, let 
(5 r denote the system with boundary conditions as in Section [2], Boltzmann weights 
B r , and rows labeled in descending order from top to bottom. Similarly, let (5 A 
denote the system with the same boundary conditions, but with Boltzmann weights 
B A and rows labeled in ascending order from top to bottom. Using this language, 
we may now summarize the results of the past two sections in a single theorem. 

Theorem 4. Given a dominant weight A for GL r+ i ; the metaplectic Whittaker func- 
tion W(t\) is expressible as either of the two partition functions Z{G T ) or Z(& A ). 

This is merely the combination of Theorems [T] and [2] together with Proposition [2J 
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5. Transfer Matrices 

Baxter considered the problem of computing partition functions for solvable lattice 
models (cf. [2]). His approach is based on the idea of using the Yang-Baxter equa- 
tion (called the "star-triangle identity" by Baxter) to prove the commutativity of 
row transfer matrices. We will show that basic properties of metaplectic Whittaker 
functions can be interpreted as commutativity of such transfer matrices, and at least 
when the metaplectic degree n = 1, the Yang-Baxter equation can be used to give 
proofs of these. 

The row transfer matrices shall now be described. Let us consider a row of vertices 
that all have the same Boltzmann weights. If B = (a\, &i, 62, c±, C2) then we use 
the assignment of Boltzmann weights in the following table. 







3 


9 


9 


0- 


3 


( 

0- 


3 

^0 


€ 


3 

fe 






e 


e 




3 




3 




3 


Boltzmann 
weight 


ax 


a 2 


61 


i>2 




c 2 



The vertical edge spins in the top and bottom boundaries will be collected into 
vectors a = (a^, ■ ■ ■ , a ) an d P = (Pn, ■ ' " > ^o)- The subscripts correspond to the 
columns which, we recall, are numbered in ascending order from right to left. For 
example, if a = (—,+,—,+,-1-,—) and = (+,— ,+,+,+, —), we would consider 
the partition function of the following one-layer system of ice: 



(22) 
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Let Vb{oc, /3) denote the partition function. Recall that we compute this as follows. 
We complete the state by assigning values to the interior edges (unlabeled in this 
figure) and sum over all such completions. Let Vb be the 2 N+1 x 2 N+1 matrix whose 
entries are all possible partition functions Vb(«,/3), where the choices of a and p 
index the rows and columns of the matrix, respectively. This is referred to as the 
transfer matrix for the one-layer system of size N with Boltzmann weights B at every 
vertex. 
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Now let us consider a two-layer system: 



(23) 
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Note that we are using two sets of Boltzmann weights B\ and B 2 for the top and 
bottom layers, respectively. We may try to express the partition function V(a,7) 
for the two layer system pictured above having top row a as in ( 122]) and bottom 
row 7 = (+, +, +, — , +, +) in terms of one-layer partition functions. However each 
one-layer system is only determined upon a choice of vertical spins lying between 
the rows. One such choice of edge spins is /3 as in the one-layer system in ([22]) . but 
we must sum over all possible choices to get the partition function of the two-layer 
system. Therefore 



which is precisely the entry V(a, 7) in the product of the two transfer matrices Vb 1 
and Vb 2 . 

Cases where the transfer matrices commute are of special interest. Indeed, this 
commutativity means that one can interchange the roles of Boltzmann weights B\ 
and B2 in ( 123]) and the value of the product of the transfer matrices is unchanged. 
Baxter considers the case where B\ = (a, a, b, b, c, c) and B 2 = (a', a', b', b', c', c') for 
arbitrary choices of a, a', b, b', c, d . However, his boundary conditions are toroidal; 
that is, the boundary edges at the left and right edges of the each row are identified 
and treated as interior edges, hence summed over. With this modification, Baxter 
proves that if A = A', where A = (a 2 + b 2 — c 2 )/2ab and A' is similarly defined with 
a', b' and c' , then the transfer matrices commute. Obtaining a sufficiently large family 
of commuting transfer matrices is a step towards evaluating the partition function, 
since by doing so one can make the eigenspaces one-dimensional. Thus the problem 
of simultaneously diagonalizing them has a unique solution and therefore becomes 
tractable. 

Let us now show how Statement B may be formulated in terms of commuting 
transfer matrices. We consider a two layer system having a layer of Gamma ice and 
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a layer of Delta ice, thus: 
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(24) 
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We use the values in f fl8|) ; in the top row, the spectral parameter is zi, and in the 
bottom row it is z^. Regarding the boundary conditions, as always the rows of ice 
must have + at the left edge and — at the right edge. Furthermore, we fix a choice 
of spins for the top edge and the bottom edge of this two-layered ice such that the 
top edge has two more — than the bottom row. In this example, the locations of the 
— along the top edge are (reading from right to left) 0, 1, 4, 6 and along the bottom 
edge, they are at 3, 4. We have labeled each vertex with • I\ and oA 2 to remind the 
reader of the Boltzmann weights that we are using. We will call this system (5 rA 
and its partition function Z(& rA ). 

On the other hand we may consider the same configuration with the roles of the 
Boltzmann weights for Gamma and Delta ice switched, as in the figure below. Note 
that the boundary conditions remain the same as in (1241) . We will refer to this system 
as 6 Ar . 




Theorem 5. Given top and bottom boundary values as vectors of spins a and 7 and 
Boltzmann weights B\ and B A as in < T7%)) . let (3 rA and <5 Ar be the systems described 
above. Then Z(& TA ) = Z(& AT ). 

We prove this by showing that the claim is equivalent to Statement B, stated as 
Theorem [3] here and proved by combinatorial means in [6] . Note in particular that 
we have reformulated Statement B as the commutativity of two transfer matrices. 

Proof. We associate two strictly decreasing vectors of integers with a and 7, which 
we call I and m. Namely, let I = (/ , ^1, h, 4 ■ • ), where the Zj's are the integers such 
that = — , arranged in descending order; m is defined similarly with regard to 7. 
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Thus, in the example (I24p above there are — spins in the 6, 4, 1, columns of the top 
row and so I = (6, 4, 1, 0), while m = (4, 3). Similarly, given any admissible state of 
the system, let /3 be the middle row of edge spins, and associate in similar fashion a 
sequence a = (ai, a 2 , ■ • ■ ) according to the location of — signs in (3. 

We observe that the sequences /, a, m interleave. This holds for the same rea- 
son that the rows of the pattern interleave in Proposition [TJ it is a consequence of 
Lemma 2 of [I]. Therefore the legal states of either system (3 are in bijection with 
the (strict) short Gelfand-Tsetlin patterns 

[4 k ■■■ 4-i It ] 

t = < di d 2 O-r-l °>r } ■ 

I rrii m 2 m r -\ J 

These are not in bijection with the terms of the sum G rA (t) appearing in Theorem [3] 
because there is no condition on the middle row sum. Rather, the states of ice give all 
possible middle row sums. However, letting Q rA (S, z) denote the Boltzmann weight 
for a state of TA ice, this may be regarded as a homogeneous polynomial in the two 
spectral parameters z\ and z 2 of our two-row system. In the notation of Proposition [21 
this monomial is zf'^'^^'z^ 1 ®~ d2 ® } where dj(t) denotes the z-th row sum in the 
short Gelfand-Tsetlin pattern above. Clearly, the middle row sum can be recovered 
from knowledge of this monomial for fixed choice of boundary conditions a and 7, 
which dictate the top and bottom row of the short pattern. A similar correspondence 
may be obtained for the Ar system whose short patterns t' are associated to the 
monomial z^ 3 ' )_dl( ' t ^zf 1 ^ '~ d2 ^ \ Of course, the boundary conditions remain constant 
whether we are using the TA or Ar system, so d (t) = d (t') and ^(t) = d 2 (i'). Thus, 
the monomials 

z *(t)-di(t), z <ii(t)-d a (t) and z *(O-*(*0 ;s <*iP0-«fe(O 
agree precisely when 

di(i) = do(t) + d 2 (t)-d 1 (fO, 

which is exactly the condition on the sum in Theorem |3J Hence we see that the 
commutativity of transfer matrices - the statement that Z(& rA ) = Z(& Ar ) - is 
an equality of two homogeneous polynomials and the matching of each monomial 
corresponds to the identity of Statement B for each possible middle row sum. □ 

6. The Yang-Baxter Equation 

The proof of Theorem [51 the commutativity of transfer matrices, uses the equiva- 
lence with Theorem |3] and hence implicitly relies on all of the combinatorial methods 
of [B] in order to obtain this result. In this section, we want to explore the extent 
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to which Baxter's methods for solving statistical lattice models, most notably the 
Yang-Baxter equation, may be used to prove the commutativity of transfer matrices. 

In our context of two-dimensional square lattice models, the Yang-Baxter equation 
may be viewed as a fundamental identity between partition functions on two very 
small pieces of ice - each having 6 boundary edges to be fixed, 3 internal edges, and 
3 vertices each with an assigned set of Boltzmann weights. 

Definition (Yang-Baxter Equation). Let R,S, and T be three collections of Boltz- 
mann weights associated to each admissible vertex. Then for every fixed combination 
of boundary conditions a, r, a, (3, p, 9, we have the following equality of partition func- 
tions: 




Recall that these partition functions are sums of Boltzmann weights over all admis- 
sible states. Hence, the left-hand side is a sum over all choices of internal edge labels 
H, v, 7, while the right-hand side is a sum over internal edge labels 4>,ip,5. Note that 
the roles of S and T are interchanged on the two sides of the equality. 

In the diagram above one vertex, labeled R, has been rotated by 45° for ease of 
drawing the systems. It should be understood in the same way as S and T - it has 
a Boltzmann weight associated to a set of admissible adjacent edge labels. However, 
vertices of this type have a distinguished role to play in the arguments that follow, 
so we use the term R-vertex to refer to any vertex rotated by 45° like R in (j25|) . 

Once equipped with the Yang-Baxter equation, the commutativity of transfer ma- 
trices, i.e. invariance of the partition function under interchange of rows, may be 
proved under certain assumptions. We illustrate the method with a three-layer sys- 
tem of ice (3 having boundary conditions and admissible vertices like those of the 
system (5 r , to give the basic idea. Suppose we wanted to analyze the effect of swap- 
ping the second and third rows in the following configuration: 
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Suppose there exists only one admissible R-vertex having positive spins on the 
right; without loss of generality we take it to have all positive spins. Then the 
partition function Z(&) for (1261) multiplied by the Boltzmann weight for the R-vertex 
with all + spins is equal to the partition function for the following configuration of 
ice. 

5 4 3 2 1 

e © © e © e 
a- ^ch 
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(By assumption, the only legal values for a and b are +, so every state of this problem 
determines a unique state of the original problem.) Now we apply the Yang-Baxter 
equation to move this R-vertex rightward, to obtain equality with the the following 
configuration. 
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Repeatedly applying the Yang-Baxter equation, we eventually obtain the configura- 
tion in which the R- vertex is moved entirely to the right. 

5 4 3 2 1 

e e e e e e 

^Q^ 
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In drawing the above picture, we have again assumed that there is just one legal 
configuration for the R- vertex having two — spins on the left, and assumed the spins 
of this R- vertex were all — . If we let & denote the system with the same boundary 
conditions as & shown in ( |26i) but with the second and third row Boltzmann weights 
interchanged, we have shown 



(29) 




Z(&) 




where Br denotes the assignment of Boltzmann weight to each configuration. In 
particular, if the two admissible R-vertices coming from the left- and right-hand 
sides of f )29|) have equal Boltzmann weights, we obtain the exact equality of the two 
configurations, i.e. the commutativity of transfer matrices. 

We now explore the possibility of obtaining a Yang-Baxter equation with S and T 
in f l25|) corresponding to the Boltzmann weights B r and B A , respectively, from f|T8|) . 
In light of our previous argument, this would give an alternate proof of Theorem [3j 
However, the Boltzmann weights in f fT8|) depend not only on spins + or — on adjacent 
edges, but also on a "charge" a mod n. Recall from Section H]that using B r weights, 
charge records the number of + signs in a row between the given vertex and the — 
boundary spin at the right-hand edge the row. Using B A weights, charge counts the 
number of — signs between the vertex and the + boundary at the left. 

In order to demonstrate a Yang-Baxter equation, we need Boltzmann weights that 
are purely local — i.e., depend only on properties of adjacent edges — so we need 
a different way of interpreting charge. We do this by labeling horizontal edges with 
both a spin and a number mod n. We declare the Boltzmann weight of these vertices 
to be unless the edge labels a, b mod n to the immediate left and right of the vertex 
reflect the way charge is counted for the given spins. For example, with B r weights, 
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a = b + 1 if the spin below a is + and a = b if the spin below a is — . Using this 
interpretation, we record the non-zero vertices for both sets of Boltzmann weights: 

m 
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The above vertices are admissible for any choice of a mod n (and the integers a + l 
are, of course, understood to be mod n as well). This means that we are generalizing 
the six-vertex model, since due to the dependence on a, each vertex has more than 
six admissible states. 

For n = 1, the charge labels on horizontal edges are trivial as the Gauss sums 
g(a) and h(a) are independent of a as evaluated in (jSJ). For this special case, it was 
shown in j4] that a Yang-Baxter equation exists with weights S and T as in (j25j) 
taken to be B T and B A from the table above. We refer the reader to [J] for the 
corresponding Boltzmann weights R for which the Yang-Baxter equation is satisfied. 
Thus we obtain an alternate proof of Theorem HJ or equivalently Theorem [3j using 
methods from lattice models. 

In general, we know from [6] that Theorem H] is true for any positive integer n. 
It would be extremely interesting to find a local relation like ([25]) similarly proving 
that the transfer matrices commute, and this is currently under investigation by the 
authors. 



7. Weyl group invariance and the Yang-Baxter equation 

Kazhdan and Patterson [131 Lemma 1.3.3] describe how the metaplectic Whittaker 
functions transform under the action of the Weyl group. This invariance — which does 
not follow directly from the description of the coefficients H given in Theorem [Tj— 
plays a key role in the proof of the metaplectic Casselman-Shalika formula for GL r+ i 
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by Chinta and Offen [9], and was the main inspiration for the Weyl group action 
in [ID]. 

In this section we restate this Weyl group invariance in terms of the partition 
functions defined in the previous sections. We content ourselves to describe how a 
simple reflection acts on the partition function. Let <7j denote the simple reflection in 
the Weyl group corresponding to the z-th simple root. We let <7j act on the spectral 
parameter z = (z 1 ,z 2 ,..., z r+1 ) by a^z) = (z u . . . , z i+1 , z h z i+2 , z r+1 ), i.e. 
the i and (i + l) st coordinates are transposed. Here the notation Z(&,z) refers 
to the partition function associated to the system &, where & is either of the two 
systems (5 r or <5 A introduced in Section @J 

Further define, for j — 0, . . . , n — 1, 

(31) pWfry) = ^V n - \lZ q q -\ yn ««i Q®(*,V) = M) f_~-i r > 

where we again use the shorthand notation of (J7J) and interpret g(0) := g(n) = —q' 1 . 
The functions P,Q are closely related to the functions t},t% of [HU Lemma 1.3.3]. 
For each 1 < i < r, we may decompose the partition function 

(32) z(e,z)= J2 

0<j<n 

where Z^\&,z) is the sum over all states 5 6 6 such that B(S,z) is equal to 
a constant times z^ 1 ■ ■ ■ z a r r +^ where a, — a i+1 = j (mod n). Then the Whittaker 
function satisfies 

(33) Z^(&, a t {z)) = P U) (z t+1 , ■ Z?\e, z) + Q«\z i+1 , z t ) • zf^G, z). 

We now consider the extent to which the functional equations ( 133]) can be inter- 
preted in the language of transfer matrices. First we consider the case n = 1. The 
decomposition on the right-hand side of ( 1321) has only one term, namely Z itself, 
since the congruence condition is automatically satisfied by all monomials for any i. 
The 2-th functional equation ( |33l) becomes 

Z(6,ai(z)) = (p( Q \z i+1 , Zi ) + Q {0 \z t+u z t ))Z(e,z), 

or 

(34) ( Zi - z i+1 /q)Z(&, z) = (-zt/q + z t+1 )Z(&, a^z)). 

Recalling the effect of crj on z defined above, the partition function on the right- 
hand side is the result of swapping the spectral parameters associated to rows % and 
% + 1 in the system &. Note that ( 134]) is not exactly the same as "commutation of 
two transfer matrices" because we do not have the identity Z(&,z) = Z(&,<Ji(z)). 
Indeed, the partition function Z is not a symmetric function, but it is very close to 
one: it is a Schur polynomial times a g-deformation of the Weyl denominator (cf. [1]). 
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Nevertheless, with assumptions as in Section [6j we may ask for a Yang-Baxter 
equation leading to a proof of (}3"1]) . That is, we seek sets of Boltzmann weights R, S, 
and T satisfying ( 125]) where S = T = B T or S = T = B A . Comparing ( 13^|) with 
( 1291) . we further require Boltzmann weights Br for the R- vertices such that 




It follows from results in [I] that we may use the following coefficients in the R-vertex 
for Gamma ice: 
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We are taking all U = — q 1 in Table 1 in [3], and observe that the order of the rows 
in this paper are opposite those in that paper. Our convention here is the same as 
in ®. ' ' 

For n > 1 the situation is more complicated, but rather suggestive. In general 
Z\ j \&) ^ Zl n ~ j) {&), so we cannot rewrite ([33]) to look like (|29l) and (jMj). However 
according to (13TT) . the denominators of P and Q appearing in the ith functional 
equation (I3UI) are equal and independent of j. For any j, they are zf — zf +1 /q. Thus, 
clearing denominators, we may rewrite (I33p as follows 

(36) pU\z i+1 , zt) ■ Z?\&, z) + qW(z i+1 , zi) ■ ^(6, *) 

where 

pW(z i+1 , Zl ) = (1 - q- l )zl +1 zr\ q U) (z l+u z t ) = g(j)(z? +l - 

Let S be a state of the system, and as before let a%, ■ • ■ ,a r+ \ be the exponents 
of zi, ■ ■ ■ , z r+ \ in B(S, z). We make the following observation. In the weights fl30|) . 
there is a contribution of Zi if and only if the charge is not augmented as we move 
across the vertex. Since (in Gamma ice) the charges at the right edge will have value 
0, it follows that the charges at the left edge will have value q where cij + Cj is the 
number of vertices in the row. Therefore 



(37) 
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and we may therefore write 

z®(e,z) = J2 B ^ s ^)- 

Cj-|-i — Ci = j mod n 

We will now explain how, with a suitable R-vertex, ( )36l) could also be interpreted 
as an identity similar to (12 9j) . but now with sets of Boltzmann weights involving 
charges. We will describe the characteristics that such an R-vertex might have. For 
simplicity, we will assume that n is odd. 

The value will depend on the spins and charges of the adjacent edges. Let us 
assume first that the spins on these four edges are all +, with charges d i+ i, di, d' i+1 , 
d\ as follows: 



di + i d i+1 




(38) d i d i 

If j = d i+ i — di and j' = d' i+1 — d[ then we require that the Boltzmann weight of this 
vertex v is zero unless j' = j or n — j mod n. Moreover in these cases we require 
that the Boltzmann weight of fl38l) is 

f p®(zi+u z i) if 3 = f mod n 
\ q^(zi+i, if j = n — j' mod n 

except when j = 0. In this case the weight will be 

p<°W.*) + 9 (0) (^+i,^) = 4 - q-'z^, 

since g(0) = —q^ 1 - 

Regarding the case where the vertex has spin — on all four adjoining edges, we 
require that the Boltzmann weight of 



di+i 




di 



is zero unless di = d i+ i = 0, in which case it is zf +1 — q~ x z™ . 

Assuming that the R-vertex has the above properties, we may now express the 
functional equation in a form similar to ( )29l) . Let us fix the vertical edge spins above 
the Zi + \ row and below the Zi row, and work with just the two relevant rows; let & 



2(> 
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denote the two-layer system consisting of just rows % + 1 and i with these boundary 
spins fixed. In order to establish (133|) . or equivalently (|36|) . it suffices to show 



p®(z iH ,z i )-Z®(6 , ,z) 



H+l 



q {j \z i+1 ,Zi) 



zt j \e',z). 



Since Z(&', z) is a homogeneous polynomial in the Zi, and since only ai and a%+i are 
allowed to vary, we have ai + aj + i equal to a constant. Since we are assuming that 
n is odd, there will be a unique pair of charges q and q+i mod n such that ( 1371) is 
satisfied, and such that q+i — q = j modulo n. 

Now let us consider the partition function of the system 






obtained by attaching the R-vertex to the left of &. From the above discussion, this 
equals 

Similarly the partition function of the system 

hl © ® 




is 



(^ 1 -?- 1 ^)^ C, ' ) (© , ^i(*)) 



The equality of these partition functions implies ( 1361) . 

At this writing, we do not know if the values of the R-vertex that we have described 
can be completed to a full R-matrix such that the appropriate Yang-Baxter equation 
is satisfied. We know that this can be done when n = 1, and since (136!) is true, it 
seems very plausible that this can be done in general. Thus we may conjecture that 
within this scheme, or some similar one, it is possible to formulate a Yang-Baxter 
equation adapted to these weights that gives a proof of (|36l) . Such a "metaplectic" 
Yang-Baxter equation might well have importance beyond the problems that we have 
discussed in this paper. 
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